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Abstract. An asymmetric exclusion process with A particles on L sites is considered 
where particles can move one or two sites per infinitesimal time-step. An exact analysis 
for A = 2 and a mean-field theory in comparison with simulations show even/odd 
oscillations in the headway distribution of particles. Oscillations become maximal 
if particles try to move as far as possible with regard to their maximum velocity and 
particle exclusion. A phase transition separates two density profiles around a generated 
perturbation that plays the role of a defect. The matrix-product ansatz is generalized 
to obtain the exact solution for finite A and L. Thermodynamically, the headway 
distribution yields the mean-field result as A^ 1 — > while it is not described generally 
by a product measure. 
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1. Introduction 

Driven-diffusive systems such as the asymmetric simple exclusion process (ASEP) have 
been extensively used to model traffic flow phenomena p.]. It is defined on a one- 
dimensional lattice in which particles can move only in one direction with respect to the 
exclusion rule which implies that at most one particle can stay at a single site. From 
a theorist's point of view these models are especially interesting in respect of phase 
transitions, such as a jamming and condensation transitions and their solvability at least 
for the steady state. Usually, the rather sophisticated models that claim to reproducing 
real traffic features are not available for exact mathematical descriptions. The more 
important are the minimalist models that lead to an understanding of the underlying 
physics. The ASEP in one dimension with periodic boundary conditions is fairly simple 
and has a uniform stationary measure [2]. There are, however, some generalizations of 
the ASEP with periodic boundary conditions that lead to phase transitions. An example 
is the ASEP with a single defect particle that can itself move forward on empty sites 
and can be overtaken by normal particles [16] . The defect can be thought of as a truck 
that moves in an environment of cars [2]. Another example is particle disorder: If any of 
the particles has an individual fixed hopping rate one might observe a phase transition 
from a fluid into a condensed phase [TJ. Finally phase transitions have been studied in 
asymmetric exclusion models in which the hopping rate depends on the empty space 
ahead. These models can often be related to the zero-range process and the interactions 
are normally long-ranged when condensation transitions appear [8j [9]. The ZRP itself 
allows for an arbitrary number of particles per site and the single-particle hopping-rate 
depends only on the occupation on the departure site [9|. This has been generalized to 
models in which more than one particle can move. A condition on the hopping rates 
has been derived for the steady state to take a simple factorized form [TO] . 

In the following section we define a simple traffic model, which is a generalization 
of the ASEP in the sense that particles can move one or two sites per infinitesimal time- 
step. We find that the system leads to oscillations in the distribution of head-ways which 
become maximal when it is impossible to move only one site if there are more empty 
sites available. Here the system evolves in special regions of the configuration space that 
gives rise to a phase transition. Although it is a very simple conserving process on a ring 
with one species of particles, no overtaking and short-range interactions it is capable to 
produce a phase transition and has a non-trivial steady state that is obtained exactly. 
We investigated also the process with parallel update and found the matrix-product 
stationary state, see [5]. 

2. The general process 

The general process we are going to investigate is defined on a one-dimensional lattice 
with L sites, enumerated I = 1,2, ... ,L. Each site / may either be occupied by one 
particle (77 = 1) or it may be empty (77 = 0). We impose periodic boundary conditions 
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and let the system evolve in continuous time. Particles can move one or two sites to the 
right according to the following rules: 

100^010, at rate pi, 
— > 001, at rate P2, 
101 011, at rate/?. (1) 

The parallel- update version of this process has been considered in [19]. Note that 
the total number of particles N is fixed due to the allowed transitions and boundary 
conditions. Some simple cases are already known: For example pi = (3 was studied in 
[23] and turned out to have a uniform stationary state. For p 2 = one finds [23J [2TJ that 
the weights (for @ > 0) are of the pair-factorized form: P(r 1; r 2 , . . . , r£) = Il£=i ^( r «) T i+i) 
with some simple two-site factors t(Ti,r i+ i). A mapping onto a mass- transport model 
shows that these are the only cases with factorized steady state (see section l3~3l) . 



2.1. Two-particle study 

In the following we consider (pQ) with only two particles on a ring. The quantity of 
interest is the un-normalized steady-state weight f(m,n), denoting that one particle is 
followed by m and the other by n holes. Analysis of the master equation shows that 
these weights obey the following second-order recursion relation 

f(m, n) = UL(pi)f(m — 1, n) + p2f(m — 2, n), for m > n and m + n = L — 2 > 3, (2) 
with the piecewise defined function 

Pi, for L even, 
1, for L odd. 



(3) 



One rate can be chosen as 1 and it is convenient to set pi = 1 to get rid of the even/odd 
dependence of the lattice size. Concluding we take for the study of the two-particle case 
u>l(Pi = 1) = 1- In terms of the functions 



l + y/l + 4p 2 



1 



the solution of 



f(m,n 



V 2 / V 2 
(for m > n and m > 1) is 

+ (2p 2 + 1 - /9)y m+n _i 



-r 



^(1-/3)^ 



l + 4p 2 



(4) 



(5) 



1 + 4p 2 

One sees that the first term only depends onn + m and therefore it is constant for 
given system size. The second term has a pre-factor (— l)™" 1 which indicates that 
in general there are oscillations. Thus the weights f(m, n) depend on the parity of 
n. The probability for a certain configuration with two particles is given by P(m, n) = 
Z m + n+2 2f{ m i n )i where for system size L the normalization Zl,2 = Dm=o /( m > L—2—m) 
is: 



PL 2(l-/3)(2p2 + l) 



l + 4p 2 



;i+4p 2 



Vl-i + 



(2p 2 + 1 - p)L 2p 2 (l - P) 



1 + 4p 2 



;i+4p 2 



yx-2-(6) 
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One might think that the feature of oscillations comes from the presence of a finite 
number of particles, so we investigate the thermodynamic limit within a mean-field 
theory. 

2.2. Mean-Field Theory 

To take into account correlations between consecutive particles we write down an 
improved mean-field theory. The quantity of interest is the probability P{m) to find a 
headway of m empty sites in front of a particle. In the context of traffic-flow models this 
is referred to as car-oriented mean- field theory (COMF) [4jJ§| The stationary equations 
read: 

(c + p 2 s)P(0) = (3P(l)+p 2 P(2), 

(c + p 2 s + P)P(1) = cP(0) + Pl P(2) + p 2 P(3), (7) 

(c + p 2 s + pi + p 2 )P{m) = p 2 sP(m — 2) + cP{m — 1) + p 1 P(m + 1) + p 2 P{m + 2), 

for m > 2, 

with the short-hand notations 

s := 1 - P(0) - P(l) and c := (3P(l) + p^s. (8) 

We now introduce the generating function 

oo 

Q(z) = £ P(m)z m . (9) 

m=0 

Summing up P(m)z m for m = . . . oo leads to a rational expression for Q(z) from which 
a singularity at z — 1 can be removed. One obtains 

= (g -ft-ttwy -^-^(Q) , (10) 

+ (C + p 2 S)Z l - (pi + p 2 )-2 - 

with iu := (px +p 2 )P(0) +p2P(l)- A useful check of this equation is Q(0) = P(0) and 
Q(l) = 1. The density p in the corresponding asymmetric exclusion process is 

00 

d z (zQ(z))\ z=1 = £(m+ l)P(m) = p'\ (11) 
This gives P(l) in terms of P(0) and p: 

on \ _ [ 2 ^(1 + P) + Pi] P(0) - (4P2 + Pl) P n9 , 

(j " (/3- Pl )(l-p)-2p 2 • U2j 

The remaining probabilities can be obtained from Q(z). The flow-density relation is 
J(p) = p{c + 2p 2 s). 

However at this stage already one equation is missing. One needs an additional 
relation between P(l) and P(0) to be able to express everything in terms of the density 
only. In fact, the missing relation can be extracted from the generating function. Writing 
the numerator of Q(z) in terms of its zeros Zq gives (f3 — p\ — p 2 )P(l)(z — Zq)(z — Zq). 

§ This formally corresponds to a mean-field theory in the corresponding mass-transport model (see 
section [373| and neglecting correlations between adjacent masses. 
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The singularity in the unit circle then has to be removed by Zq or Zq for Q(z) to be 
analytic [4J. This leads to the missing relation between P(l) and P(0) [T5] . 

We restrict ourselves here to the case where the coefficient of z 2 in the numerator of 
Q(z) vanishes, since there the mean-field predictions can be written in a very compact 
form and the main-features that we want to display are contained. Consider (3 = pi+p 2 '- 
a headway's total rate of change is independent of its length. The condition is weaker 
than the condition for a factorized state (every configuration is equally probable if 
Pi = (3). Substituting this into ( fTUl) and demanding that the denominator has the same 
zero as the numerator gives the missing relation: 

p 2 P(l) 2 = Pl P(0) (P(0)[1 - P(0)] - P(l)) . (13) 

For Q(z) one gets the simple expression 

= pPm 2 (14) 

pA 2 — pAz — z 2 

with A := P(0)/P(1) and p := P1/P2, which can nicely be expanded to obtain P(m). 
Figure [1] shows the mean-field distribution for two different choices of parameters 
in comparison with computer simulations. While P(n) decays rapidly, one sees the 
characteristic even/odd oscillations which are well reproduced by mean field. 




Figure 1. Headway distribution P(n) from a computer simulation with L = 1000 in 
comparison with mean field for p2 — 1. Left: p\ — 0.1, (3 = 1.1 and p = 0.1. Right: 
pi = 0.2, (3 = 1.2 and p = 1/3. 



2.2.1. The Fibonacci case Remarkable is the case (3 = p\ + p 2 = 2 with pi = p% = 1. 
Then the probabilities P(m) are given by the Fibonacci numbers: 

P(m) = P(0) (^j F m+1 . (15) 

Here one obtains from <$j5§: P(1)/P(0) = (^5 - 4P(0) - l)/2 and relating P(0) to the 
density gives 
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Figure 2. Headway distribution P(n) from mean-field for the Fibonacci case p\ = 
P2 = 1 and (3 = 2 for p — 1/8 and 1/2 compared with computer simulations. 



Figure [2] shows the headway distribution for the Fibonacci case. 



2.2.2. The choice p± = In (jTUl) one sees that also -P(l) = reduces the numerators 
degree to one. It has a zero at z = —pziPi + £>2) _1 - We are interested in p 2 > 0, so 
take pi — 1 without loss generality. Demanding that the denominator has the same 
zero yields (as the only physical solution) p\ = 0. The generating function reduces to 
P(0) (1 + (1 - P(0))z 2 + (1 - P(0))V + . . .)). So P(l) = has the consequence that 
P(2n + 1) vanishes generally and the process is realized only on the even sublattice in 
mean field. The relation to the density is 

P(0) = T^—y (17) 
l + p 

and the flow simply reads 

J(p) = 2p[l- P(0)] = 2P(1 ~ P) ■ (18) 

l+p 

These results are completely equivalent to the usual ASEP where now particles always 
move two sites and the density is appropriately rescaled. 



3. Maximal Oscillations: The choice p\ = 

We consider the process 

100^001, at rate 1, 

101 -v 011, at rate/?, (19) 

Physically, this is the case where every particle tries to move as far as possible with 
regard to his maximum velocity. This is the limit in which the oscillations in the form 
of even/odd effects become maximal. This type of process evolves into special regions 
of the configuration space. We try to find the full solution for finite iV and L from the 
matrix-product ansatz. 
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3.1. Unified solution for finite number of particles and sites 

The state of the periodic system can be expressed by the sequence of headway: 
{ni, . . . , riN }. We look for a product state of the form: 

N 

F(n 1 ,n 2 , ■ ■ ■ ,n N ) = Tr JJ Q nil , (20) 

Ai=l 

where Q Hi is an operator representing particle /i followed by n M holes. It turns out that 
this ansatz (with a certain trace- like operation Tr=tr| |.| | to be specified below) yields the 
correct steady state, provided that the involved operators fulfill the following quadratic 
algebra: 

02i01 = 02i+l, (21) 

02i+l02j+2 ~ 02i+l02j = 0, (22) 
02i02j+2 — 02%02j = P02i+2j+2, (23) 

02i02j+3 — 02i02j+l = P02i+2j+3, (24) 

02i+i02j+i =0, forz,j>0. (25) 

We just note here that this can be proven by use of the canceling-mechanism [15]. The 
corresponding tagged operators read: 

=0o- /31, (26) 

02(i+l) = 02{i+l) + 02i, (27) 

01 = 0u (28) 

02i+3 = 02i+3 + 02i+l- (29) 



It is important to emphasize that (1251) implies that the system can not support more 
than one odd gap. This can be understood directly from the dynamical rules ffT9]) . 
The fact that the transition 100 — > 010 is forbidden (pi = 0) has the consequence that 
the number of odd gaps decreases with time: A configuration C(. . . l[any odd number 
of 0s] 101 1 . . .) moves with conditional probability (3 into a configuration with two odd- 
valued gaps less (creation of even gaps), while odd gaps can not emerge. These processes 
appear until there remain either no more odd gaps (L — N even) or exactly one odd gap 
(L — iV odd). In the latter case this means physically that the probability for odd gaps 
is of order 1/N and thus tends thermodynamically to zero as predicted by mean-field. 

One has to be careful with a unified description for an arbitrary number of odd 
gaps in the system. The operators are mathematical objects that can in our case be 
written as matrices whose components are itself matrices (compare [T2l [5] ) . Thus it is 
not obvious how to generalize the trace operation. The straight-forward generalization 
to a sum of the traces of the matrices on the main diagonal can here not be applied: 
This trace operation and therewith the weight for certain configurations can incorrectly 
give zero. The problem of the trace operation for periodic systems has previously be 
pointed out [6]. The matrix relation fT25|) implies that (02i+i) 2 = for all i. For a unified 
description the operators 02i+i had to be non-vanishing nilpotent matrices. However 
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for this equation to hold all the eigenvalues of G 2 i+i must be equal to zero. Therefore 
also tr £?2«+i = which implies for example through relation (j2ip that also tr Q 2 M\ = 
which is untrue. Now this argument can be used successively to see that for every iV 
the straight-forward trace operation can not be applied. First, the algebra (I21H25I) can 
be simplified: 

g 2i+x = £ l A and Q 2l = 0V (30) 

with new operators S and T>. Then (12111251) reduces to 

VS =V + S, (31) 

AS = A, (32) 

(3VA = A, (33) 

A 2 = 0. (34) 

Introduce the two by two matrices 1 = |1)(1| + |2)(2|, 23 = |1)(2| and let further be 
£ = E®1, V = D®1, and A = A <g> 22, with matrices E, D, and A. Then one 
can interpret the trace operation in (120]) as F{n^ n 2 , ■ ■ ■ , n N ) = tr| 1 iT^Li ^njl> with the 
help of the matrix norm ||M|| = maxjjmjj to obtain always the correct weights. A 
different method is for example a parity-dependent matrix representation [15]. However 
beyond the question of how to obtain the correct matrix element it is most convenient 
to consider both cases separately. 

3.2. Expectation values and phase transition 

For even number of holes only even gaps occur in the stationary state. With the particles 
always making two steps (100 — > 001) the stationary process is completely equivalent 
to the ASEP. All configurations with even-length gaps have the same stationary weight 
(the matrices T> and 8 are sufficient to describe the steady state and can be chosen as 
numbers). For the normalization we find 

L{L-M-\)\ S 

Z L,N - N\M\ ° L - N > 2M - \ 6b > 

This is easily interpreted combinatorially: For the first particle one has L possible ways 
to place it on the lattice. One can then think of distributing N — 1 particles and M 
hole pairs into N + M — 1 boxes to obtain the above expression. The flow (from site i) 
is the expectation value 



J = Ml - r l+1 )(l - n+2 )) + (3{ri{l - T l+1 ) n+2 )) = (n) - ( n r m ). (36) 



With (n) = f and (7-^+1) = N/L ■ (N - 1)/(L - M - 1) this yields asymptotically the 
mean- field current (|T8l) . The exact form of the velocity for finite system size reads 

l-p(2 2 n \ 

^ 2 TT7A 1 + TTn + - + TlTW + -)- (37) 



Headway oscillations and phase transitions for diffusing particles 



9 



For odd number of holes exactly one odd gap occurs in the steady state. If one considers 
only stationary configurations the relation (1341) becomes redundant. The non- vanishing 
steady-state weights from ( 120]) can simply be written as 

"iV-l 

F(2ni,2n 2) ...,2n JV + l) = tr JJ E n »D E nN A. (38) 

Here we referred to the particle with the odd gap in front as particle N. Due to the 
translational invariance this can be done without loss of generality. The underlying 
algebra reduces to 

DE = D + E, AE = A, and (3D A = A. (39) 

This is the algebra for the ASEP with a single defect particle [H] for the case a — 1. 
This process is defined by transitions 10 — > 01 at rate 1, 20 — > 02 at rate a = 1 and 
12 — > 21 at rate (3. In fact, the stationary states of both processes are completely 
equivalent. In our process even gaps of length 2n are mapped onto gaps of length n in 
the defect ASEP (00 becomes 0). Now consider the single odd gap with the particle to 
its right. It can be written in the form 00 00 ... 01. Again 00 is mapped onto and 01 
becomes the defect 2. One can check that under this mapping indeed ([19]) recovers the 
transitions of the defect ASEP. Note that a = 1 takes care of the (physically reasonable) 
fact that in our process particles with an even or odd gap to their left move forward at 
the same rate. 

As for the defect ASEP the partition function can be calculated: 
LfN + M\™ (N + M- 1\ (l -/?\ m_1 . 

From this expression correlation functions can be derived as above. It turns out that 
the flow is related to the normalization by 

J = 2 - z l-W +(3 L - NZ L-W-\ (41) 
L — 2 Zl,n L — 1 Zl,n 

We have calculated the finite size expansion for the velocity in the case (3 = 1 and 
obtained 

1-p ( 5/2 (1 + 3 n+1 )/4 \ 

The important thing is that the correction is of order 1/(L + N) which of course holds 
also for (3 ^ 1. 

To summarize in both cases (even and odd number of holes) the velocity of particles is 
given by 

• = 2 Tr| + °(rb)- (43) 

Just the special form of the correction differs for even and odd number of holes. Consider 
now explicitly a finite number of particles. The weights can be obtained easily from the 
algebraic rules. For two particles the weights are of the form f(2l + 1,2m) = 1 + 1(3. 
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For three particles we find f(2l + 1, 2m, In) = 1 + (I + n)/3 + + In) (3 2 , and so 

on. From these formulae the probabilities P(n) to find a gap of n sites (n = 0, 1, ... 4) 
can easily be calculated. The results are shown in figure [31 One sees the remarkable 
change of the distribution by adding a single empty site to the system. Note that the 
probability for an odd gap is of order 1/N and thus independent of L. Therefore the 
qualitative form of holds also for L — > oo. 




Figure 3. Headway distribution P(n) for p2 = 1, pi = 0, and (3 = 1 for N = 3 and 4 
and L = 102 and 103. 



For the thermodynamic limit, one might directly make use of the density profile [21 [16] . 

It is interesting that the single 'excess hole' created by the dynamics leads to effects 
also in the thermodynamic limit. Although the probability P(n) as predicted by mean 
field is thermodynamically exact, the system reaches no product-measure steady state, 
since around the defect a nontrivial density profile is formed. Beyond that a phase 
transition takes place equivalent to a transition in the defect ASEP [TH]. In terms of the 
different density (p in the (defect) ASEP becomes here 2p/(l + p)) the critical density 
is p c = (3/(2 — (3). In the following let us in analogy refer to the 01-pair as the defect. 
Since we have a = 1 the phase diagram of the defect ASEP reduces to a single line. 

• For p > p c the defect behaves as the other particles. In front of the defect the den- 
sity profile decreases exponentially to its bulk value. The density behind is constant. 

• For p < p c the defect is similar to a second-class particle p2] that lowers the average 
speed of the other particles. The density profile decays algebraically to the bulk 
value. Behind the defect the density is decreased and the profile increases in the 
same way to its bulk value as in front. The profile is the limit of a shock profile 
with equal densities to the left and right. 

From the relation to the defect ASEP one can obtain the probabilities P(2n + 1) for 
odd headway. For example P(l) is related to the probability p_ in [16] to find a particle 
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directly behind the defect. Since in our process the defect can be any of the N particles 
one has 



[p(mp) 



4p 2 



(5N(l + pf 
2p 



for p < p c , 
for p > p c . 



(44) 



[N(i + P y 

Figure H] shows -P(l) scaled with N versus the density for (3 = 2/3, so that the phase 
transition happens at p c = 1/2. Depicted are the analytic formulae from (T441) together 
with a computer simulation for L = 1000 with N increased in steps of AiV = 25. 




Figure 4. -P(l) versus p for (3 — 2/3. See text for details. 



3.3. Relation with generalized zero-range processes 

In the previous section we expressed the steady state of the system through the set of 
gaps between particles. This formally corresponds to a mapping onto a model in which 
a site occupied by particle /i becomes site /i and the gap n^_i to the left becomes the 
'mass' on site fi. The ansatz fTSOj) then becomes site-oriented. The corresponding process 
comprises N sites and M := L — N particles. In the process obtained from fTTOT) by this 
mapping one or two particles may leave a certain site with rates: 

!1, for I = 2,m > 1, 
R f 7 1 1 (45) 

p, for / = 1, m = 1. 

However this is a special case of the class of generalized zero-range processes introduced 
by Evans, Majumdar and Zia [10J and referred to as 'mass-transport models' (MTM). 
They derived a necessary and sufficient condition for the steady state to factorize, i.e. 
P(mi, ni2, ■ ■ ■ , mjy) oc n„=i f(m u ). The condition on the chipping functions 7 reads 

nam: 

7(/M = ro( ^ (m -" .^ ' = 1.-,™ («) 
f(m) 
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where w(l) is an arbitrary non- negative function of I. For the process (14 5 p the situation 
is slightly more special. Here we have a factorized steady state only for even mass and 
thus it is not predicted by ( )46|) . In our rather singular case, involving vanishing single- 
site weights, (H6|) is not fulfilled, as is easily seen for 7(1 |m). To understand this, let us 
consider a slightly more general case. Assume general single-site weights that vanish for 
odd masses as it is the case for fj45l) . Considering the master equation shows that the 
restrictions for the chipping function to recognize the factorized state for M even, are: 

free, for m + I odd, 



j{l\m + l) 



0. for /, m both odd, 



x(l)f(m) 



for L m both even. 



f(m + iy 

However the choice for f(m) in turn implies that the total mass M = J2i=i m % is even. 
Otherwise the normalization 

z NM = E<M M -E^')n/K) ( 4 8) 

{mi} \ j / i=l 

would vanish. The solution can be used to obtain new solvable models: The process 
fT45|) leads to a matrix-product state for the choice of M for which the system can not 
reach the factorized state. This suggests that the general model ( 1471) with odd particle 
number (defined trough the weights for even M) may also lead to a matrix-product 
state by special choice of the free parameter in ( 1471) . Of course these arguments can 
be generalized to other choices of vanishing weights [15J. Beyond that, the single-site 
mass-distribution should equal thermodynamically the result obtained from the case 
where it is factorizable, since in the infinite system a local perturbation changes the 
density profile but not the single-site distributions. This way one can obtain the exact 
distributions also for cases where the steady state has not generally a product measure. 
Some focus has recently be attended to two-species zero-range processes and conditions 
for a factorized steady state [9]. We consider model parameters that violate this 
condition. Denote the number of first-class particles on site I by n\ and the number 
of second-class particles as m ; . The rates u(n, m) and v(m, n) at which first and second- 
class particles move are taken as 

u(n, m) = 1, for n > 1, (49) 

v(m, n) = (3 [1 — 0(n)] , for m > 1. (50) 

The motivation for this choice is simply that the steady state of the MTM (and 
equivalently (|T9l) ) corresponds to the case of a single 2-particle here (pairs of particles 
in the mass-transport model are mapped onto 1-particles in the ZRP, and the single 
excess particle is mapped onto the 2-particle). The resulting algebra then also becomes 
a consequence of ( 1391) . 
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4. Conclusion and outlook 

To summarize, a simple traffic model that generalizes the ASEP with periodic 
boundaries has been considered. In this continuous-time process, particles can move one 
or two sites to the right. This mimics a larger maximum velocity of cars that is typical 
for discrete traffic models [3]. Some preliminary results have already been published [20J. 
We considered choices of hop rates that are not appropriate for modeling traffic but are of 
theoretical interest. A mean- field theory in comparison with simulations showed that the 
headway distribution of particles show even/odd oscillations. We considered a special 
choice of the parameters (see ( fl9i) ) for which the oscillations become maximal. The 
matrix-product ansatz was generalized to obtain the exact solution for finite number N 
of particles and sites L. Thermo dynamically a crucial phase transition appears: Particles 
move at their desired maximum velocity v max which leads to the formation of headway 
being a multiple of v max . If there remains smaller headway, then it moves through the 
system against the direction of motion as a sort of defect. This defect (being somewhere 
in the system) changes locally the density profile but not the headway distribution of 
a single particle. Therefore the system is not described by a product measure although 
the asymptotic headway distribution is given by the mean- field result. Instead the 
solution is the matrix-product state for the ASEP with a single defect [16J. Note that 
the process with parallel dynamics leads to a similar stationary state with additional 
attraction between particles and hole-pairs and can also be solved [5]. 
We established a relation between non-ergodic exclusion processes with higher velocities, 
generalized zero-range processes (ZRP) and defect systems. This way one can calculate 
exact quantities without knowing the exact density profile. For future work it is 
interesting to investigate the connection between systems with creation and annihilation 
of 'defects' and generalized ZRP to be able to handle ergodic dynamics without parity 
dependence. 

The condition for a simple factorized state in the ASEP considered here is that the rate 
at which a particle moves a certain number of sites is independent on the headway. This 
condition holds also for parallel dynamics [22] and higher v max . We investigated a slightly 
more general parameter line, where the rate at which a particle can change its headway 
does not depend on its length which allows for oscillations in the headway distribution. 
The mean-field assumption leads to a remarkable agreement with computer simulations. 
We pointed out the Fibonacci case whose corresponding mass-transport model (on two 
sites) has a factorized steady state. The single-site weights become Fibonacci numbers 
s(n) = F n+ \. The recursion relation ([2]) becomes f(m,n) = f(m — l,n) + f(m — 2,n), 
for m > 2, which can be expressed as a matrix-product state f(m, n) = tr (DE m DE n ) 
with DEE = DE + D. However neither the factorization nor the matrix recursion 
hold for more than two particles. So this case remains unsolved. A comparable good 
agreement with mean field has been observed previously for the ASEP with shuffled 
update [24J where also the two-particle state is factorizable. It would be interesting to 
find out whether the mean-field result in these models is generally exact although the 
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process has not a product measure. Beyond that, the connection between the solvability 
for two particles and N particles is still an open problem. Here one might be able to 
profit from knowledge in equilibrium statistical mechanics |25j . 
We would like to thank A. Schadschneider for helpful discussions. 
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